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Abstract 



We study the slow decay of the steady-state autocorrelation function C{t) 
in a stochastic model of deposition and evaporation of trimers on a line 
with infinitely many conservation laws and sectors. Simulations show that 
C{t) decays as different powers of t, or as exp(— t^^^), depending on the 
sector. We explain this diversity by relating the problem to diffusion of hard 
core particles with conserved spin labels. The model embodies a matrix 
generalization of the Kardar-Parisi-Zhang model of interface roughening. In 
the sector which includes the empty line, the dynamical exponent z is 2.55 ± 
0.15. 
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It is well known that conservation laws strongly affect the decay in time 
of fluctuations in equilibrium. For example, in magnetic systems, the rate 
of decay of the spin autocorrelation function is quite different depending 
on whether or not the spin dynamics conserves magnetization [|l|. What 
happens with more than one - indeed an infinity of - conservation laws? We 
address this question in this paper by studying the decay of autocorrelation 
functions in the steady state of a simple stochastic model of deposition and 
evaporation p|, ||. This model has been shown to possess an infinite number 
of independent constants of motion [Q. Here we show that it exhibits a very 
rich variety of slow relaxations: depending on the initial conditions, we can 
get a very large number of power-law decays, or even stretched exponential 
decay. We present numerical evidence from Monte Carlo simulations for these 
different kinds of decay, and also provide an analytical understanding of this 
remarkable diversity in terms of the diffusive motion of hard core particles 
with spin. Our model can be viewed as a generalisation of the well-known 
Kardar-Parisi- Zhang (KPZ) problem ^ of roughening of a one dimensional 
interface, where the scalar-height variables in the KPZ model are replaced by 
a matrix- valued variables (say 2x2 complex matrices). In the steady state 
reached after starting from an empty lattice, we find that the density-density 
autocorrelation function decays as a power law in time, and the dynamical 
exponent z ~ 2.5, very different from the KPZ value z = 3/2. The model 
thus includes a new universality class for nonlinear one-dimensional evolution 
models. 

In the deposition-evaporation (DE) model under consideration, there is 
a site variable n^, taking values or 1, at each site i of a line, 1 < i < L; 
it may be thought of as an occupation number variable of a lattice gas. We 
define pseudospin variables Si = 2ni — 1. The time evolution is Markovian, 
specified by the following rates: In a small time dt, a triplet of spins at 
adjacent sites i,i + l,i + 2 can flip simultaneously only if Si = Si+i = 5*^+2- 
The rate is e if the spins were originally -1, and e' if they were originally 
+ 1. These rates satisfy the detailed balance condition corresponding to a 
non-interacting lattice gas Hamiltonian with chemical potential |£n(e/e'). 
However, the long-time steady states in the present model are nontrivial, as 
the dynamics is strongly non-ergodic. We have shown that the total phase 
space of 2^ configurations breaks up into a large number Nl of mutually 
disconnected sectors, and Nl increases as [{V5 + l)/2]^ for large L. 

This decomposition of phase space is a consequence of the existence of an 
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infinity of conservation laws. A compact representation of these is provided 
by the "irreducible string" (IS), defined as follows: Regard the configuration 
as a string of up and down spins. Scan the string left to right, and stop at 
the first triplet of parallel spins encountered. If no such triplet is found the 
string is the IS. If a triplet is encountered, simply delete it, reducing length of 
string by 3, and repeat the process. The IS which finally results is a constant 
of motion, and different sectors of phase space correspond to different IS's 

i- 

Under time evolution, an initial configuration evolves into the steady state 
of the corresponding sector. For instance, if e = e', in the steady state all 
configurations having the same IS as that of the initial configuration occur 
with equal probability. It is then interesting to ask how a dynamical quantity 
such as the steady state autocorrelation function 

Ci{t) = {ni{t + to)ni{to)) - {rii^ (1) 

varies from one sector to another. Figure 1 shows the behaviour of C{t) 
obtained from Monte Carlo simulation studies in a number of different rep- 
resentative sectors. The data depicted in the figure pertain to four sectors in 
which the IS constitutes a finite fraction of all sites, and one sector in which 
the fraction is zero. These sectors correspond to (1) a random initial config- 
uration (2) an IS formed by repeating [t||] L/6 times (3) an IS formed by 
repeating [fl] L/4 times (4) an IS formed by repeating [ttitii] L/12 times 
(5) the null sector (vanishing IS). Lattice sizes L = 120,000 were used and 
averages over 30 — 60 different histories were taken. From the figure, we note 
the following points: (i) if the lattice is divided into 3 sublattices A,B,C, the 
autocorrelation function Ci(t) depends, in general (e.g in Sectors 2 and 4), on 
the sublattice to which site i belongs, (ii) in several sectors and sublattices 
Ci(t) decays as a power law ~ t^^; the power 9 is sector dependent. In Sector 
1, 6* ~ 0.25; in Sector 2 and on two sublattices of Sector 4, 6^ ~ 0.5; in Sector 
5, 6' ~ 0.6. (iii) in some sectors (e.g. Sector 3 and one sublattice in Sector 4), 
the decay is faster than a power law, and fits well to a stretched exponential 
form. Evidently, there is a very wide range of possible behaviours. Under- 
standing the source and nature of this dynamical diversity is one of the main 
points of this Letter. 

Consider time evolution in a sector labelled by an IS with elements = 
{ai, a2, ■ ■ ■ , ae}, where the length £ is a nonzero fraction p = i/L of the total 
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length L of the lattice. Let C{t) be the full string corresponding to the 
configuration at time t. If we apply the deletion algorithm to C(t), some 
characters are eventually deleted, others not. If Xj{t) is the location of the 
j**^ site (counting from the left) that remains undeleted, we may look upon 
the set {Xj{t)} as the positions of C. interacting random walkers on a line; 
DE dynamics induces the time evolution oi{Xj{t)}. The walkers cannot cross 
each other > Xj{t), for all j, all t) and they carry a conserved spin 

label Sxj{t) = «j for all t. Under time evolution Xj{t) changes in jumps, by 
a multiple of 3 spaces at a time. 

In the steady state, the joint probability distribution Prob{{Xj}) is pro- 
portional to the number of different configurations of the lattice, consistent 
with the positions of walkers being {Xj}. The substring between sites Xjj^i 
and Xj should be a string reducible to a null string, and all such configura- 
tions are easily enumerated. The result is 

Prob{{X,}) = AT n g{X,+, - X,) (2) 

3=0 

where is a normalization constant, and Xq = 0, X^+i = L + 1. The 
separation probability g{r) can be computed using the generating function 
method of 0]. The result is g{r) ~ r~^/^exp(— Kr) for large r, where k is the 
reduced pressure, which depends on on the density £/L, and tends to zero 
as iS-jVf for small ijL. We see that {Xj\ constitute the slow modes of the 
system as they evolve diffusively, and are linked to conserved quantities. If 
the typical relaxation time of a string of length A reducible to the null string 
varies as A^, then for times t ^ {L/iy we can assume that all the degrees 
of freedom other than {Xj} are in instantaneous local equilibrium. Thus 
g{r) ~ exp{—V{r))^ where the effective interaction V{r) between walkers is 
attractive, and increases linearly with separation if k 7^ 0, and as In r if 

K = 0. 

While the logarithmic part of the interaction is crucial for the under- 
standing of sectors in which the walker density £/L — 0, it appears not to be 
very important if i/L, and hence k, is finite. In the latter case, we are led to 
consider a simpler system: i random walkers on a line of length L, with each 
walker carrying a spin aj which is unchanged under dynamics. Each walker 
jumps left or right by 3 steps only if no other walkers intervene. The point 
is that this simpler system of hard core random walkers with spin (HCRWS) 
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has the exactly the same conservation laws as the original deposition evap- 
oration model. Thus we are led to conjecture that the long-time behaviour 
of Ci{t) in a particular sector of the DE model is essentially the same as the 
spin-spin autocorrelation function Di{t) in the HCRWS problem with the 
corresponding spin sequence. 

It is then easy to understand the sublattice dependence of Cj(t) observed 
in simulations. As each element a„ of the IS moves by multiples of 3 lattice 
spacings, it stays on the same sublattice, say A. The long-time behaviour 
of Ci(t) ioT i & A is governed only by the subset {a^'j^i on sublattice A, 
independent of {a^'ls and {an'}c, despite the fact that the elementary step 
of deposition or evaporation couples all sublattices strongly. 

For the HCRWS model, mean squared fluctuations in the number of par- 
ticles between points and r grow as t^^^ in the steady state 0. To test 
for a similar behaviour in the DE system, we monitored cT^(t) = {[N{r,t) — 
A^(r, 0)]^) where N{r,t) is the length of irreducible string between points 
and r at time t. This quantity is roughly equal to number of walkers to the 
left of r, but remains well defined even in the sector where there are no walk- 
ers i = 0. Figure 2 shows Monte Carlo data for (j'^{t) in a number of sectors. 
For all sectors with i/L finite, cr^(t) is seen to grow as t^^^, lending strong 
support to the conjecture of equivalence of the DE and HCRWS models. In 
the sector where the IS vanishes, on the other hand, o"^ is found to grow as 
t"^^ with f3 ~ 0.19. If the irreducible string is very short, then k is very small, 
and fluctuations in the separation (Xj+i — Xj) between near neighbours in 
the diffusing, interacting lattice g large. The low value of the exponent 
f3 implies that these fluctuations are slowly decaying. 

In the HCRWS model, the spin-spin autocorrelation function in the steady 
state is defined as 

D,it) = {S,{t + to)SM)-{S,)^ (3) 

where Si(t) is spin at site i, and an unoccupied site is considered to have zero 
spin. Di(t) can be written in the form 

+ 00 

D,it)= Yl G{m\tH{m) (4) 

m=—OD 

Here, given that in the steady state site i occupied at some time Iq, G{m\t) 
is the conditional probability that it is occupied at time [t + to), and that the 
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difference in the particle labels at these times is 3m. This term is independent 
of the spin configuration {ctj} of the IS. The second term 7i(m) is the average 
value of akak+3m, averaged over different values of k in the IS consistent 
with the condition that particle k can occupy site i (i.e. k — i mod 3). The 
different dependences of Di(t) on time in different sectors comes entirely from 
the different dependences of 7,j(m) on m. 

Using the known equivalence [7,8] of long-time HCRW dynamics to a 
stochastic harmonic model, G(m|i) for large t is given by 

G(m\t) ^ — ^e-"^'/252 (5a) 

where 

Using eq. (4) and (5) we can determine the behaviour of Di{t) in the HCRWS 
model in different sectors. 

Sector 1 (random initial condition): In this case ji{m) is significant only for 
very small values of m. For m fixed, G ~ t~^/^ for large t, whence Di{t) also 
varies as Note that this result is also true for other sectors in which 

the correlations in the IS decay fast with distance. 

Sector 2 [tii]: In this case ct^ = ^fc+sm for all m. Hence Di{t) reduces 
to a density-density correlation which decreases as t^^^"^ for large t. The 
same behavior is expected in all periodic IS where the magnitization in each 
sublattice is nonzero. 

Sector 3 [jj,]: In this case 7i(m) = (—1)™", and this implies that Di{t) has 
the stretched exponential form exp[— (t/r)^/^] at large time at all sites. The 
same behaviour would occur whenever the IS is periodic with zero net mag- 
netization on each sublattice. 

Sector 4 [TtiT-li]- sublattices A and C we have = ak+sm, and thus as 
in sector 2, we have Di(t) ~ On sublattice B, 'yi{m) is (—1)™, and as 

in sector 3 we get stretched exponential decay. 

If the sector is such that correlations in the IS imply that %{q), the 
Fourier transform of 7i(m), varies as \q\'^ for small q, then Di{t) would vary 
as 

It is evident that the HCRWS predictions for Di{t) accord very well with 
the Monte Carlo results for Ci{t) shown in Fig. 1. The occurrence of stretched 
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exponentials in the [||] sector, in particular, explains why earlier attempts 
[0] to fit finite size data to power laws yielded anomalously low values of the 
dynamical exponent z. 

For any configuration {Si\ of spins on the line, we define the matrix 
variable at site i by the equation 

h{t) = \lA{SM (6) 
i=i 

where A{1) and A(— 1) are 2x2 complex matrices given by 

where uj = exp(27ri/3) and x is a real parameter. Then as = A{—lY = 

1, it follows that /^(t) is a constant of motion |Q. The evolution of the 
variables {/^(t)} is Markovian, governed by local transition rates. If we start 
with an empty line, Ii{t = 0) take simple values, and get progressively dis- 
ordered as time increases. For large x, the ratio In Tr Ir{t)/ln x behaves in 
qualitatively the same way as N{r,t). At t = 0, the variance cr^ of N{r,t) is 
zero, and grows to a value proportional to r in the steady state. If N{r, t) 
obeyed a stochastic evolution equation of the Kardar-Parisi-Zhang type P, 
this variance would grow as t^/^. Our numerical results (Fig. 2) clearly 
rule out this possibility. From Figs. 1 and 2 we see that the decay of the 
autocorrelation function is characterized by ~ 0.6, while correlations in 
the length of irreducible string up to a given point grow anomalously slowly, 
with (3 ~ 0.19. Static correlation functions in this sector are characterized by 
power law decays as well. For instance, if we define an indecomposable sub- 
string as one which can be completely reduced, but which cannot be written 
as the concatenation of smaller substrings reducible to the null string , the 
probability of occurrence of an indecomposable substring of length r varies 
as r~^/^ for large r. This is consistent with a random-walk picture of fiuctu- 
ations in nearest-neighbour distance, and thus a static correlation exponent 
X = 1/2. The dynamical exponent z = x/ P thus ~ 2.6. A numerical 
diagonalization study @| of the spectrum of the relaxation operator for finite 
systems with L < 30 yields a value 2.55 ± .15. This value of z is nonstan- 
dard, indicating that the dynamics in this sector is characterized by a new 
universality class. In this sense, our model is different then earlier attempt 
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1^ to generalize KPZ by allowing the height function to be an component 
vector; it was argued that z = 3/2 independent of in c? = 1. 

To summarize, we have shown that the deposition-evaporation model 
exhibits great diversity in the manner in which the autocorrelation function 
decays. There is strong evidence for the conjecture that in sectors with finite 
IS density i/L, the dynamical behaviour is essentially the same as that of 
the spin correlation function of a system of hard core random walkers with 
the corresponding spin sequence. The diversity of relaxational behaviours in 
the DE model is linked to differences in correlations in spin patterns in the 
IS which labels that sector. In the sector composed of completely reducible 
configurations, the dynamics is characterized by a critical exponent z ~ 2.5, 
indicative of a new universality class. 

MB acknowledges the hospitality of the Condensed Matter Group at the 
International Center for Theoretical Physics, Trieste where part of this work 
was done. 
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Figure Captions 



Figure 1: Variation of the decay of the autocorrelation function Ci{t) in 
Sectors 1-5 (defined in text). Sector 1 (open triangles); Sector 2, sublattice 
A (+), sublattices B and C (filled triangles); Sector 3 (filled circles); Sector 
4, sublattices A and C (open squares), sublattice B (open circles); Sector 5 
(filled squares). The straight lines show power law decays with 9 — 0.25 
(top), 0.50 (middle), 0.62 (bottom); the curve shows a stretched exponential 
(~ exp(— t^/^)) decay. 

Figure 2: Time-dependence of the mean squared fluctuations in the length 
of irreducible string N{r, t) up to point r. The four sets of data (from top) 
pertain to Sectors 1-4 (shifted vertically for clarity) while the fifth pertains 
to Sector 5. The straight hues have slopes 0.50 and 0.38. 
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